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LIMIT THEOREMS FOR RANDOM WALKS 


ALEXANDER BENDIKOV, WOJCIECH CYGAN, AND BARTOSZ TROJAN 


Abstract. We consider a random walk St which is obtained from the simple 
random walk 5 by a discrete time version of Bochner’s subordination. We 
prove that under certain conditions on the subordinator r appropriately scaled 
random walk St converges in the Skorohod space to the symmetric a-stable 
process B°‘. We also prove asymptotic formula for the transition function of 
St similar to the Polya’s asymptotic formula for 


1. Introduction 

It is well known and due to Polya [15] that the transition density pa(x, t), 0 < a < 
2, of the one-dimensional symmetric a-stable process has the following asymptotic 
expression ^ 

(1.1) Pa{x,t) ~ , as t\x\~^ —)■ 0 

where 

( 1 . 2 ) 

For multidimensional symmetric a-stable process similar asymptotic formula has 
been found by Blumenthal and Getoor [7]. Both results were obtained by using 
the method of characteristic functions and the scaling property of a-stable distri¬ 
butions. An elegant proof of the Polya-Blumenthal-Getoor asymptotic formula 
was proposed by Bendikov [2] . Bendikov used the fact that the symmetric a-stable 
process can be obtained from the standard Brownian motion by means of the sub¬ 
ordination, the notion introduced in the theory of Markov semigroups by Bochner 
[ 8 ], 

From the probabilistic point of view, the a-stable process = {Bf : t > 0), 
0 < a < 2, is obtained from the Brownian motion B = [Bt : t > 0) by setting 
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= Serf where the subordinator a°‘ = (erf : s > 0) is independent of B non¬ 
decreasing Levy process such that 

E(e-^<) = 6-"^“''", A>0. 

From the analytical point of view, the transition density Pa{x, t) of B°‘ is obtained 
as a time average of the transition density p(x, t) of B, i.e. 


Pa{x,t) 



p{x,s) di/t“(s) 


where is the distribution function of the random variable erf. In particular, the 
minus infinitesimal generator of the process B°‘ is (—A)“/^ where A is the classical 
Laplace operator. 

In this article we study a certain class of random walks obtained from finite 
range random walks, in particular the simple random walk, by means of discrete 
time subordination recently developed in [4] . We prove that any random walk from 
our class, being appropriately scaled, converges in the Skorohod space to the a- 
stable process 5“. We also prove that a discrete space and time counterpart of the 
Polya’s asymptotic formula (I.l) holds true. 

To illustrate the main results of the paper we consider S to be the simple random 
walk in a discrete space and time counterpart of the Brownian motion in 
Let 0 < a < 2 and S°‘ be the subordinate random walk defined via the functional 
calculus of (minus) generators of discrete time Markov semigroups 




Here P and Pa are the transition operators of S and S°‘, respectively. In other 
words. S'" coincides in distribution with St-^ where the discrete subordinator r“ is 
uniquely determined by the formula 

/ / \ ^ 

E(e-^^") = n - ^ A>0. 

In Sections 4 and 5 we prove the following two theorems (or even more general 
versions of them). 


Theorem A. For each 0 < a < 2, the sequence of random elements 

(n-i/"S[)),] : f > 0) 

converges in the Skorohod space to B°‘. 

Let {ei,..., ej} be the standard basis for K'’* and || • II 2 the Euclidean norm. 


Theorem B. Let 0 < a <2 and Pa{x,n) be the transition function of S°‘, then as 
n and ||a ;||2 tend to infinity we have 
(i) j/n||a:||^“ tends to zero then 


1 


d 

ipc‘ix+ 

j=i 


ej,n) + 2pa{x,n) + Pa{x 


- ej,n)) 
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(ii) ifn\\x\\ 


—a 

2 


tends to infinity then 

Pa{x, n) ^ PaiO, n) ~ 

The constants in Theorem B are given explicitly by the formulas 

'd + a\ . /'Tra\ 


Cd,. = a2“/^7r-'^/^-id-“/^r(|)r(^) sin (^) 


and 


7^,^ = (2^)d/2£(l±^, 

’ ^ ^ r(i + d/2) 

Let us observe that Cd^a has the factor d~°‘^'^ which is not present in (1.2). This is 
a consequence of the asymptotic behaviour of the simple random walk on 

The main ingredients in our proofs are the continuity principle of superposition in 
the Skorohod space due to Whitt [21], and the asymptotic formulas for the discrete 
time subordinators. In Section 2.3 we prove the following theorem. 


Theorem C. Let Fn{t) = P(t„ < t). Whenever nt “ tends to zero, 


1 - Fn{t) 


r(l-a/2)- 


Let us mention that Theorem C can be regarded as a version of the global 
large deviation result for random walks with increments belonging to the domain 
of attraction of a stable law, see e.g. [9, 14] and references therein. 

In this work we consider more general subordinators. Namely, for any Bernstein 
function ^ which varies regularly at zero we construct related subordinator r and 
random walk St, and we establish r-versions of Theorems A, B and C. 

Theorem A (as well as the more general version in Theorem 4.1) was first ob¬ 
tained by A. Mimica in [13]. Furthermore, he also proved the converse of the 
theorem and treats the case a = 2. While Mimica applied the method of charac¬ 
teristic functions, our approach is different and uses superposition in the Skorohod 
space. 


2. Preliminaries 

2.1. Finite range random walks. Let S = {Sn : n G N) be a random walk driven 
by a probability measure p and let p^^'> = p{- ,n) be its n-fold convolution. We say 
that {S,p) G ii the following properties hold 

(i) p is supported by a finite set V C 

(ii) for each x G there is n G N such that p{x, n) > 0; 

(iii) E«eZ‘iP(^)-^ = 0- 

Generally, a random walk from the class is not aperiodic. Thus, we let r to be 
its period, that is the greatest common divisor of all n such that p{0, n) > 0. Then 
the space decomposes into r disjoint classes i?o, •. ■, Rr-i where 

Rj = {x G : p{x,j + kr) > 0 for some fc > O}. 

For instance, if S is the simple random walk starting at the origin, we have r = 2 
and i?o is the class of points which S visits in even number of steps. 
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Before we state the asymptotic formula for p(x, n) we need to introduce a qua¬ 
dratic form on defined by 

(2.1) {Qu,u) = ^ p{v){v,uf. 

We note that for the simple random walk Q = dT^I where / is the identity matrix. 
Let Q~^ be the inverse of Q. We also fix a norm on R'^, ||a:||^ = {Q~^x,x). By [18, 
Corollary 3.2 and Remark 1], we obtain that for every j G {0,..., r — 1}, n G N 
and X G Rj such that ||x|| < if n = j (mod r) then 

(2.2) p{x,n) = (2Trn)~i {detQ)~^ -|-(!l(n“^)-I-(!l(n“^||x||^)^ 
otherwise p(x,n) = 0. In particular, we obtain 

Corollary 2.1. For every e > 0, and x gU^, 

r — 1 

^p(x, n + i) = (27rn)“^ (det Q)~^ -I- 0[n~^) + C>(n“^||a:||^)^ 
j=o 

provided ||a;|| < . 

Let us comment about the asymptotic of p(x,n). For aperiodic random walks 
with mean zero to get (2.2) one may use [12, Theorem 2.3.11]. A more refined 
asymptotic is recently developed by Trojan in [18]. It is valid in a larger region and 
applicable to random walks with non-zero mean and not necessary aperiodic. 


2.2. Discrete time subordination. For continuous time Markov processes, sub¬ 
ordination is a useful procedure of obtaining new Markov processes. The later may 
differ very much from the original one but its properties can be understood in terms 
of the old process. For example, in this way the symmetric stable process can be 
obtained from the Brownian motion (see e.g. [2]). 

From a probabilistic point of view, a new process (Yt : t > 0) is obtained from 
the process (Xt : t > 0) by setting Yt = X^t where the subordinator (as : s > 0) is a 
non-decreasing Levy process taking values in (0, oo) and independent of (Xt '■ t > 0) 
(see e.g. [10, Section X.7]). From an analytical point of view, the transition function 
ha(x, t) of the new process is obtained as a time average of the transition function 
h(x,B,t) of the original one, i.e. 


hcr(x,B,t) 



h(x, B, s) diyt(s) 


where vt is the distribution function of the random variable at. Subordination was 
introduced by Bochner in the context of semigroup theory (see [10, footnote, p. 
347]). Formally, the minus inhnitesimal generator B of the process (Ft : t > 0) is a 
function of the minus infinitesimal generator A of the process (Xt : t > 0), that is, 
B = tjj(A) for some function il). 

A discrete time version of subordination, where the functional calculus equation 
B = ^/’(A) serves as the defining starting point, has been considered by Bendikov and 
Saloff-Coste in [4]. To describe the procedure, let R be the operator of convolution 
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by p, a probability density on Then L = I — R, where I is the identity operator, 
can be considered as minus the Markov generator of the associated random walk. 
For a proper function ip we may define a subordinate random walk S'^ as the random 
walk process with the generator —ip^L). As it has been shown in [4], the appropriate 
class consists of Bernstein functions. 

Let us recall that a smooth and non-negative function ip : [0,oo) —> [0,oo) is 
called a Bernstein function if for all n G N and a; > 0 

< 0. 


For a Bernstein function ip there are constants a,b > 0 and a measure v on [0,oo) 
satisfying 


/ min{l, s} di/(s) < oo, 

[0,oo) 

ip{x)=a + bx+ ( (1 — dz/(s). 

J [ 0 , 00 ) 


and such that 
(2.3) 

Given a Bernstein function ip with ip{0) =0 and ■0(1) = 1 for fc G N we set 

b if /c = 1, 

0 otherwise. 

Then the random walk S'^ with the generator —ip{L) has the probability density 


(2.4) 


^) = T\ i + 

a[0,oo) 


(2.5) 


Pip{x) = '^p{x,k)c{ip,k), 
k>l 


see [4, Proposition 2.3]. Formula (2.5) has the following probabilistic interpretation. 
Let {Rk : fc G N) be a sequence of independent and identically distributed integer¬ 
valued random variables also independent of S and such that P(i?i = k) = c(ip^ k). 
Then for Tn = Ri + ■ ■ ■ + Rn we have 

(2.6) P(t„ = fc) = ^ c{ip,ki) ■ ■ ■ c{ip,kn). 

ki + ...+kn=k 

In particular, the law of is the same as the law of Sr„ and is given by 

(2.7) p^(a;,n) = ^p(a;, fc)P(r„ = fc). 

k^n 


2.3. Asymptotics related to subordinators. For a given {S,p) G =5^ and a 
Bernstein function ip, the relations (2.6) and (2.7) defining the transition function 
p^{x,n) of the subordinate random walk S'^ are rather complicated. In the course 
of study, we are going to restrict ourselves by considering Bernstein functions ip 
regularly varying at zero of index 0 < q;/2 < 1, such that '0(0) = 0 and '0(1) = 1. 

Let us recall that a function / defined in an interval (0,a), a > 0, is regularly 
varying of index P at zero if f{x) = x^£{l/x) where £ satisfies 
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for each A > 0. The function £ is called slowly varying at infinity. For the detailed 
exposition of the theory of regular variation we refer the reader to [6]. Here, we 
only mention Potter bounds (see [16], see also [6, Theorem 1.5.6]). If £ varies slowly 
at infinity then for every e > 0 and C > 1 there exists xq > 0 such that for all 
x,y>xo 

(2.8) £(y) <C£{x) ma.x{y/x, x/yY- 

The main result of this section is the following theorem. 


Theorem 2.2. Set Fn{t) = P(r„ < t). If n and t both tend to infinity in such a 
way that n'ip{t~Y tends to zero then 

nip{t~Y 


1 - Fn{t) 


T{l-a/2y 


Proof. The statement of the theorem may be expressed as follows: for any sequence 
{tn '■ n € N) tending to infinity and satisfying 


(2.9) 
one has 

( 2 . 10 ) 


lim 


lim nV'(£„ = 0, 

n—>-oo ^ ^ 

l-F^tn) 1 


n-ipitn^) r(l — a/2) 

Let us observe that the Laplace transform of the function Fi is equal to 

pCO 

(2.11) £{Fi}(A)=/ Fi(y)e-^^dy = A-i(l-V'(l-e-^)). 

Hence, for any n gN 

(2.12) £{1 - £„}(A) = A-i(l - (1 - y{l - e-^))”). 

Next, we define 

J^n{x) = [ 1 - Fn{y) dy, 

Jo 

for X > 0. Then, by (2.12) 

F„(x) = r dFYv) < e l\-yl^dTYv) 

JO JO 

< eC{dFn}{I/x) 

(2.13) <Cxml}{l-e-^/^). 

In what follows C denotes an absolute positive constant which may vary form line 
to line. 

For a given sequence (£„ : n G N) satisfying (2.9) we define a sequence of tem¬ 
pered distribution on [0,oo) by setting 


A„(/)= 


1 


^ Jo 


f{x)J^n{tnX) dx 
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for any / G 5([0, oo)). Let us recall that the space 5([0,oo)) consists of Schwartz 
functions on R restricted to [0, oo) and 5'([0, oo)) consists of tempered distributions 
supported by [0,oo), see [19] for details. 

We claim that there is C > 0 such that for any / G 5([0, oo)) and n G N 

(2.14) |A„(/)| < CsupKl +a;)^|/(a;)|}. 

x>0 


Indeed, by (2.13) we have 

^ POO 

|A„(/)| < 2;|/(x)|V'(l - dx. 

To bound the right-hand side, let us recall that ipiX) = A“/^£(l/A) for some 0 < 
a <2. By (2.8), there are TV > 1 and C > 0 such that for n > TV 

Thus, we can estimate 

^- 1/2 t “^/2 

[ ” x\f{xM{l - dT < Csup {|/(x)|}tl“/ 2 +i )/2 / " ^ drr 

'^(in ) Jo x>0 Jo 

< C'sup{|/(x)|}. 

T>0 


Let Xn = -1/ {tn ■ log (l - 1/Tn)) and 

An{x) = T„(l - 

For tn < X < Xn we have 

1 < An { x ) < 

Again, by (2.8), there is TV > 1 such that for n > TV 

i{An{x)~'^tn) < £(Tn)A„(a:)(^““'^^)/^, 

thus 

'ip{An{x) ■ / N(a/2-l)/2 < ^ 

An(x)^(tn^) - 

Since xAn(x) < 1 whenever x > 0, we obtain 

xlf ( x ) jT / j ( An ( x ) ■ t-^} dx 


< 


*-1/2 


x|/(x)|A„(x) dx 


< Csup{(l -hx)2|/(x)|}. 

ai>0 


Finally, if x > x^ then Aji(x) < 1. Therefore, by monotonicity we get 


i/>(A„(x) < ^(t^^). 
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Hence, 


^ pOO PCX) 

J -t-^) dx < J a;|/(a;)|dx 

< CsupKl +a;)^|/(x)|}, 

x>0 

which finishes the proof of the claim (2.14). 

Let us observe, that (2.14) implies that the family of distributions (A„ : n S N) is 
equicontinuous. Next, we calculate A„(/t-) for fr{x) = r > 0. The integration 
by parts yields 

1 ^oo ^ 

(2.15) Anifr) = - , , _is / e ™J'„(t„x) dx = - £{dJ'„}(r/t„), 

ntnW[tn Jo - 


We need the following observation. 


nTtn'>l’{tn^) 


Claim 1. Let (a-n : n € N) and {bn : n G N) be any two sequences of positive 
numbers both tending to zero and such that 


Then 


lim -— = 1. 

n—Kx> bn 


,■ V' On . 
TTTT = 

n^oo 1p{bn) 


Indeed, by (2.8), for every A > 1 and e > 0 there is (5 > 0 such that for all n G N, 
if 0 < On, < (5 then 

V’(a„) < Ai/’(5„)max{(a„/6„)^+", (a^/fen)^”"}. 

Hence, 


'lP{bn 


< A, 


and since A > 1 was arbitrary, the limit is bounded by 1. Analogously, we show 
the reverse inequality. 

Now, since 

1 - e-^Ar. 


lim 

n—>00 'T/t,, 


by Claim 1, we get 


^(1 


_ p-T/tr, 


lim -—s— 

n-yoo ^{tn ) 


= 1 , 


= 


Hence, by (2.12) and (2.15) 

1 r°° 

(2.16) lim Anifr) = ^-rr / e-^^x^-°‘/^ da;. 

^ ' n^oo r(2-a/2)7o 

Finally, density of the linear span T of the set {fn ■ t > 0} and equicontinouity of 
(A„ : n G N) allows us to extend the limit in (2.16) over all of / G 5([0, 00 )) giving 

liin An{f) = — - j— / /(a;)a:^““/^ da;. 

n->oo r (2 - a/2) Jo 
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For the completeness of our presentation we give a sketch of the proof that T is 
dense in 5([0, oo)) (see e.g. [19, 20]). By the Hahn-Banach theorem it is enough 
to show that if A € 5'([0,oo)) and A{(j>) = 0 for all cf> G T, then A is the zero 
functional. Assume that A vanishes on T and let CA{z) = A(e“^ ), > 0, be 

the Laplace transform of A. Since A = 0 on T we get £A(A) = 0 for A > 0. But 
the Laplace transform is analytic in the half-plane ^z > 0, whence CA{z) = 0, for 
^z > 0. Using the connection between Laplace and Fourier transforms, 

A(^) = lim^ £A(A-f J.^), in 5'([0, oo)), 

we obtain that A = 0, i.e. A is the zero functional as desired. 

Next, we claim that 


(2.17) 


lim 


^ n,(^n) 


ntni’{tn^) F(2-q;/2) 

For a fixed e > 0, we choose (j)+ G 5([0, oo)) such that 0 <((>+< 1 and 


(j)+{x) = 


1 for 0 < a; < 1, 
0 for 1 -I- e < T. 


Then we have 


C^n r^n pOO 

{tn) = / dj'„(s) < / 4>+{s/tn) dJ^n{s) < / (()+(s/t„) d J"™ (s). 

^0 Jo Jo 


Hence, by the integration by parts we get 

Yl (t 




% < -A„(<(.;). 


Therefore, we may estimate 
. -1 


lim sup 


< 


pCO 

■ / s^““'^^(()()_(s) ds 

Jo 


ntn4’{tn^) r(2 - a/2) Jq 

1 1 

= rTi-s"“/^(()+(s) ds < —-- j— 

r(l-a/2)7o r(2-a/2) 


(l + e)^-“/2. 


Similarly, by taking </_ G 5([0,oo)), 0 < </_ < 1 satisfying 


cl)_{x) = 


1 for 0 < x < 1 — e, 
0 for 1 < x. 


one can show that 


lim inf ■ 


J~n (jn ) 


> 


ntntp{tn^) r(2-a/2) 

Since e was arbitrary, we conclude (2.17). 


(1-e) 


l-a/2 
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To prove (2.10) we follow the line of reasoning from [6, Theorem 1.7.2]. Let 
e > 0. The function 1 — Fn{s) is non-increasing therefore 


(2.18) 

Fnit) - Fn{{l - e)t) = 

/■ (1 




and 



(2.19) 

Fn{{l + e)t) - Fnit) = J 


By (2. 

17) and regular variation of i/. 



Fn{s)) ds > et{l 


Fn{s)) ds < et{l 


Fu{t)) 


Fn{t)). 


lim 


Fn{{l - e)tn) 


(1 - 

r(2-Q;/2) ■ 


Hence, (2.18) implies 


lim sup ^-4^ < 


1 l-(l-e)i-“/2 


n^oo" r(2 — q;/2) 

Similarly, by (2.19) we get 


,. . , 1 - F„{tr,) ^ 

limmi-^—7^— > 


1 (1-he)i-“/2 _ 1 


n->oo nip(tn^) T{2 — a/2) 
Finally, by taking e tending to zero we obtain (2.10). 


□ 


3. The domain of attraction 

Let {Xn : n G N) be a sequence of independent identically distributed random 
variables with common distribution /r. The distribution /i belongs to the domain of 
attraction of the a-stable law if for every n G N there are a„ G and G (0, oo) 
such that the random variable 

n 

( ^3 ~ ^n)/bn 

converges in law to a-stable random variable. As in the Section 2.3 we assume that 
Ip is the Bernstein function regularly varying at zero of index a/2, 0 < a < 2, such 
that ip{0) = 0 and V'(l) = 1- The main result of this subsection is the following 
theorem. 

Theorem 3.1. The law of Sf belongs to the domain of attraction of the symmetric 
a-stable law. 

We start with two elementary lemmas. Let $ and denote characteristic 
functions of p and , respectively. 

Lemma 3.2. For all f G 

$^(e) = l-i/(l-$(e)). 
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Proof. Notice that may take complex values, thus we need to use the holo- 
morphic extension of the Bernstein function ip. For any z G C such that Jiz > 0 we 
may write (see [17, Proposition 3.5]) 


tp{z)=a + bz+ ( (1 — e 

J [0,oo) 

The function il}{z) is continuous for IRz > 0 and holomorphic for > 0. Thus, for 
any z G C such that < 1, we have 




z) = l-b{l-z)- f (1 - du{t) 

J fO.oo') 


= bz 


' [0,oo) 


' [0,oo) 
_ fn 


^>l 


= bz + ^ [ e ‘t" diy(t) 

n>l 


n>l 


where the coefficients c{ifj,n) were defined by (2.4). On the other hand 

xgz ** k>i 

= y]c(V',fc)(4>(o)^ 

k>l 

which finishes the proof. □ 

Lemma 3.3. For u > 0 and 5Rz > 0, 

|l/’(u(l + z)) — tp{u)\ < \z\'ljj{u). 

In particular, we have 

i/>(u(l + ie)) = '!/’(u)(l + 0(e)). 


Proof. From the representation of the Bernstein function we get 

+ z)) - iIj{u) = buz + f e““‘(l - dz/(t). 

[0,oo) 

Since |1 — e“^| < jzj for Jiz > 0, we obtain 


|'0('«(1 + z)) —'0(u)| < 6u|z| + jzj f ute dz/(t) 

J [0,oo) 

< 6u|z| + jzj f (l - e“"‘) di'{t) = |z|'0(m), 
J [0,oo) 

where we used the inequality e^ >1 + x. 


□ 



12 


ALEXANDER BENDIKOV, WOJCIECH CYGAN, AND BARTOSZ TROJAN 


Proof of Theorem 3.1. Let a„ = 0. For (6„ : n G N) we choose a sequence of real 
numbers such that 

(3.1) lim = 1. 

n—¥oo ' ' 

We are going to show that for every f G 

lim =exp{- 2 -“/ 2 (Q^,^)“/ 2 } 

n—^■oo ^ ' 

where Q is defined in (2.1). Since for some C > 0 and for all a: G M 


1 — cos X - — 


< 


then for Q G [—tt, tt)"^ we can estimate 

1 
2 


5i(l - $(6»)) - i(Q6»,6») < ^piv) 1 - cos{v,e) - ^{v,d}^ 

vGV 


(3.2) 

Next, since 


<Ci||0||(Q0,0). 


'^p{v) ■v = 0, 

vGV 


and there is C > 0 such that for all x G R 

|sinx — x| < C\x\^, 

we get 

p{v) (sin(x, 9) — {v, 0)) I < E I sin(u, 0) — (u, 0) | 


vGV 


(3.3) 


vev 

<C2\\0\\{Q0.0). 


Therefore, by Lemma 3.3 and estimates (3.2) and (3.3), we obtain 


|i/>(l - $(0)) - '0(K(1 - «>(6')))| < 


3$(6») 


i/>(5R(l-$(0))) 


5R(l-$(6l)) 
<4C2||0||i/>(5R(1-$(0))). 


(3.4) 

provided ||0|| < (4Ci)“^. Now, let us fix ^ G R'^. We have 


where 






lim = 1, 

n—^oo 


Since 
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there is TV > 0 such that for n > N, > 1/2. Hence, by Lemma 3.2 and 

(3.4), we get ^ 

Log$^(^/6„) = Log(l - 'ip{l - 4’(C/&n))) 

= -1/(1 - <i>(^/6„)) (l + 0{\ij{l - $(e/6„)) d) 

= _^(SR(1 _ $(f/&„)) (l + 0(6-1) ^ _ ci>(e/6„))))). 


Next, by (3.2), we have 


thus, by applying Claim 1 


— (QC,e)/(2fo2) 

we obtain 


1 , 


lim 

n—foo 


i/(jR(l - $(g/6„))) 


lim 

n—¥cci 


tp{bn'^) 


= 2-“/2(Qe,0“^'- 


Finally, by (3.1), we conclude 

lim nLog4>^(e/&n) = -2-“/2(ge,0“/"- 


□ 


4. Functional Limit Theorem 

Without loss of generality we prove the functional limit theorem in the one¬ 
dimensional case. For simplicity we consider random walk {S,p) G such that 
Sn = X]r=i where Xi are independent and identically distributed random vari¬ 
ables taking values in Z and satisfying 

(4.1) EXi = 0 and Var(Xi) = 1. 

We consider the discrete subordinator (r„ : n £ N) defined by equations (2.4) and 
(2.6) and we set tq = 0. We assume that "0 is a Bernstein function regularly varying 
at zero of index a/2, 0 < a < 2, such that 1/(0) = 0 and i/(l) = 1. Let 5''^ be the 
subordinate random walk and let B be the Brownian motion and B°‘ the symmetric 
a-stable process in In this section we aim to prove the following theorem. 

Theorem 4.1. Assume that (6„ : n £ N) is a sequence of positive numbers such 
that 

lim 711/(6“^) = 1. 

n—^oo 

Then the sequence of random elements '■t>0) converges in the Skoro- 

hod space to the random element B°‘. 

Let us recall that a sequence : n £ N) of random elements converges to 

a random element X in the Skorohod space D([0, oo),]R) equipped with the Ji- 
topology if the following two conditions hold (see [11, Theorem 16.10 and Theorem 
16.11]): 

(i) The finite dimensional distributions of converge to the distribution 
of X. 


^Log denotes the principal value of complex logarithm. 
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(ii) For any bounded sequence {Tn : n G N) of X^^^-stopping times, and any 
sequence {hn : n G N) of positive numbers converging to zero we have 


lim 

n-^oo V 


-X, 


(") I 


Trt \ — 


> e = 0. 


Our approach to prove Theorem 4.1 is based on the continuity of a composi¬ 
tion mapping in the Skorohod space. More precisely, if a sequence of random ele¬ 
ments converges to {X, r) in the product Skorohod space D([0, oo), R) x 

D([0 , oo), [0, oo)), and if the trajectories of the limiting process Xt are continuous, 
and the trajectories of the process Tt are non-decreasing then the composition of 
the processes 

M 

converges to the process Yt = in D([0, oo), R), see [21, Theorem 3.1] for details. 
We also use the following invariance principle, for any S & satisfying (4.1) we 
have 


(4.2) 


lim = Bt 

n->oo y/n 


in the Skorohod space D([0,oo),R), see [5, Theorem 16.1]. In view of the above, to 
prove Theorem 4.1 we show the convergence of appropriately scaled subordinator 
(r„ : n G N) and then we combine it with (4.2). 

We start with the following elementary lemma. 

Lemma 4.2. Let (6„ : n G N) be a sequence of positive numbers sueh that 

lim ntfibf^) = 1. 

Then (6“^r„ . n G hi) converqes vn distribution to (0^2) stable law. 

Proof. Let F’i(t) = P(ti < t). We observe that the argument given in the proof of 
[10, Theorem 2(b), Section XIII.6], shows that if 

(4.3) l-Xi(t) ^ 


r(a/2) 


t-°‘/'^e{t), 


as t tends to infinity, then b^^Tn converges in distribution to (a/2)-stable law. 
Hence, we just need to verify (4.3). Since by (2.11) 

£{l-i^i}(A)~A“/2-i£(l/A), 

as A tends to zero, the asymptotic (4.3) is a consequence of the Tauberian theo¬ 
rem together with the monotone density theorem (see e.g. [6, Theorem 1.7.1 and 
Theorem 1.7.2]). □ 

For 0 < ,5 < 1, let us denote by cr^ = (erf : t > 0) the one-sided /3-stable 
subordinator which is uniquely determined by the following relation 




= e 
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Theorem 4.3. Let (bn : n G N) be a sequence of positive numbers such that 

(4.4) lim ntjj{b~^) = 1. 

n—^oo ^ ' 

Then the sequence of random elements t > O) converges in the Skorohod 

space D([0,c»), [0, oo)) to the random element . 


Proof. Let First, we show that the finite dimensional distributions 

of converge to the finite dimensional distribution of Since processes 

have independent increments it is enough to prove convergence of their one 
dimensional distributions. Let us fix t > 0. We write 



7 * U 

On 0[nt] 


Since cr“^^ has (a/2)-stable law, by Lemma 4.2 we obtain that converges 

in distribution to From the other side, by the uniform convergence theorem 

for regularly varying functions we get 

lim 

n—)-oo bn 

Hence, : n G N) converges in distribution to Because t°"l'^a^^'^ and 

cx /2 

(T(' have the same distribution the conclusion follows. 

Next, let (T„ : n G N) be any bounded sequence of ^("'^-stopping times, and 
{h„ : n G N) any sequence of positive numbers converging to zero. We are going to 
show that for any e > 0 

(4.5) lim f( IX^T^h - I > e) = 0. 


Since 

^T„ + h„ ^T„ 


X 


(n) y (") _ 7,-1 1 

T„+/i„ W[nT„+nIi„] ''■[nT„]j7 

he same distribution as b~^T[nhn]+jn where 
jn = [nTn + nhn] - [nT„] - [nh„]. 


We notice that jn G {0,1}, thus 

- Pnh„]) > ^) = > Kf) < IP(ti > 6„e), 

hence 

Ji^P(6-i(r[„,,„]+,„ - > e) = 0. 

We conclude that to prove (4.5) it suffices to show that 

(4.6) lim F{b-'^T[nhJi > e) = 0. 

n —' L j / 


Indeed, suppose the sequence (nh„ : n G N) is bounded by m G N. Then 
lim P(5“ > e) < lim P(rm > bnc) = 0. 

n—¥cci L j / n—¥cci 
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Otherwise, let {rij : j S N) be any increasing sequence, and {rij^ : k G N) its 
subsequence with the property that : k G N) is strictly increasing. We 

set 

Observe that by (4.4) and a regular variation of ip we have 

lim mkip{€~^b-^ ) < lim nj^iP{e-^bpi^ ) = 0. 

fe—>00 k ^ o :) 

Therefore, we can apply Theorem 2.2 to estimate 



Hence, 

lim F{b~^T[„h^] > e) = 0, 

n^oo 11/ 

which implies (4.5). □ 

Proof of Theorem 4-T Since S and r are independent, we apply [5, Theorem 3.2] 
and combine it with (4.2), and Theorem 4.3 to get that 

lim {bn^S[b^t],b~^T[nt]) = 

n—voo L 1 / \ / 

in the product space D([0,oo),M) x D([0,oo), [0,oo)). Moreover, the trajectories 
of B are continuous whereas the paths of are increasing, whence we can apply 
[21, Theorem 3.1] to get convergence of the related composition, 

lim b-hf. = 

n^oQ 

in D([0, 00 ), M). Since B a /2 has the same distribution as 5“, the proof is hnished. 

□ 

Remark 1. In paper [13] author considered the converse statement in Theorem 

4.1, i.e. that convergence in the Skorohod space implies regular variation of the 
Bernstein function related to the discrete subordinator. Moreover, the case a = 2 
is treated therein as well. 

5. Asymptotics of transition functions 

Let {S,p) be a random walk from the class . We assume that ip is the Bernstein 
function regularly varying at zero of index a/2, 0 < a < 2, such that ip{0) = 0 
and ip(X) = 1. By S'^ we denote the ^/-subordinate random walk defined in Section 

2.2. Let p^[x,n) be its transition function. In this section we are going to show 
the asymptotic behaviour of p^{x,n) as ni/(|]x|]“^) tends to zero or infinity while 
n and ||x|| tend to infinity. 
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5.1. The strong ratio limit theorem. We start by proving the asymptotic of 

PV.( 0 ,n). 

Theorem 5.1. Let p be the one-step probability for random walk S in the class 5^. 
Then 

as n tends to infinity, where 

(detQ)-i/ 2 . 

Proof. By the Fourier inversion formula and Lemma 3.2 we can write 

(5.1) p^(0,n)=^^^ (l-i/>(l-$( 6 l)))”d 6 l 

where Sid = [— 7 r, 7 r)'^. We claim that 
Claim 2. For z € C, \z\ < 1, 

|1 — '!/'(! — z)\ = 1 if and only if z = 1. 

Indeed, as in the proof of Lemma 3.2, for z G C, \z\ < 1 we have 
l-'0(l-z) = '^c{il),k)z^, 

k>l 

where the series converges absolutely for | 2 :| < 1. Since 

|l - i/i(l - z)\ < '^c{tp,k)\z\'^ < '^c{ip,k) = 1 , 

fc>l k>l 

the equation \\ —— z)\ = 1 implies \z\ = 1. Now, if |zo| = 1 and satisfies 
|1 — ' 0(1 — zo)! = 1 then there is f G M such that 

y]c( 0 ,fe)z^ = e*‘, 

k>l 

thus 

^k _ it 

Zq — e 

for all /c > 1 with c(0, k) > 0. In particular, zq = Zq, thus zq = 1- 
Next, by (2.8), there is xq > 0 such that for all 0 < a;, y < xq 

(5.2) 0(x) < 20(y)max{(x/y)“'^'^,(x/y)^“/'‘}. 

Take e > 0 satisfying 

where Ci and C 2 are constants from (3.2) and (3.3). Since |1 — 0(1 — $(0))| = 1 if 
and only if d G 27rZ'^, there is 0 < 77 < 1 such that 

|1 - 0(1 - $( 6 »))| < 1-77 
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for all 6 G = {9 G ■ ||0|| > e}, which implies 

{hy L d0 < (1 - vT. 

Next, we consider the integral over the e-ball centred at the origin. Let (a„ : n € N) 
be a sequence dehned by 



By the change of variables we can write 


a 


— d 
n 



(l-'0(l-$(6»)))"d6» 


||<e/a„ 


(l-i^(l-<i>(a„0))”de 


We are going to calculate the limit of the integrand for a fixed ^ G By (3.2) 
and (3.3), we have 

1 - d>(a„e) 1 

™ aW^,0 2’ 

thus. Claim 1 implies 


lim n^(l — ^(anO) 

n—^oo ^ ' 


lini ~ ^(Qn^)) V’(Qn(Q^.0/2) 


because n ^ Therefore, 

lim nLog(l-V'(l-‘l>(a„0)) 

n—¥cio 

Using polar coordinates, we can calculate 

I“-( - -if = 

thus to finish the proof, we need to show that the integrand (l — i/)(l — $(a„f)))” 
on G K'^ : ll^ll < e/a„} has an integrable majorant. Indeed, by (3.4) and the 
choice of e, 

|1 _ ^(1 _ ci,(a„^))| < 1 _ _ ci>(a„^))). 

From the other side, by (5.2) followed by (3.2), we get 

i^(5ft(l - <i>(a„0)) > Cn-i min {(Q^, e)“/^ {Q^, 
for some constant C > 0. Therefore, if ||^|| < e/a„ then 

- d>(a„^))P < exp C'min{(Q^,^)“/^,(g^,03“/^}y 
Since the left-hand side is an integrable function on K'^, the proof is completed. □ 


The following Corollary provides a variant of the strong ratio limit theorem. 
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Corollary 5.2. Let p be the one-step probability for random walk S in the class 
. Then 


p^{x,n) ~ p^(0,n) 

as n'!/;(||a:||“^) tends to infinity. 


Proof. By the Fourier inversion formula and Lemma 3.2, we may write 

\p^{x,n) - p^{0,n)\ < |l --0(1 - - l| dd 

<C||a;||/ |1-1/^(1-$(0))|”||0||d0. 

JSSd. 

Now, essentially the same argument as in the proof of Theorem 5.1, shows that 
there is C > 0 such that for all n G N and a; G 

\p^{x,n) -p^(0,n)| < C\\x\\'tl)~^ {n~^)~, 
which, together with Theorem 5.1, implies 

|p^(a;,n) -^^-(0,11)1 < Cp^{t),n)\\x\\'if~^ . 

To conclude the proof we need to show that \\x\\‘^'tp~^ (n~^) approaches zero when 
n'0(||a:||“^) tends to infinity. To see this, let us observe that is regularly varying 
of index 2/a at zero. By (2.8), there is d > 0 such that if n~^ < 6 and '0(||a;||“^) < (5 
then 

thus 

\\xffj-^{n-^) < 2{nf;{\\x\\-^))~^^‘", 

which finishes the proof. □ 


5.2. Asymptotic of p.^ as n'0(||x|| tends to zero. 

Theorem 5.3. Let p be the one-step probability for random walk S in the class 5^ 
having a period r. As n'0(||a;||“^) tends to zero we have ^ 

(5.3) (do + H- *p^{x,n) - rCd,an\\x\\~‘^'tjj{\\x\\~'^), 

where ||x|p = {Q~^x,x), and 

Ct.=o2V^,--/2-i(dete)-‘/A(|)r(l^) si„(^). 

Proof. First, we observe that by Tonelli’s theorem 

OO OO 

p^^^ *p^{x,n) = *p{x, fc)P(T„ = k) = ^p(x, k + j)P(T„ = k), 

k=l k=l 

thus 

CO r — 1 

{So +p + H- *p.ip(x,n) = ^^p(a;, fc + j)P(T„ = k). 

k^l j^O 


^^0 is the Dirac mass at the origin. 
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Given ko > r, we consider a sum 

ko r — l 

Si{x,n) = ^^p{x,k + j)F{Tn = k). 
k—1j—0 

Using the bounds of Alexopoulos [1, Theorem 1.8], 

p{x, k) < Ck~'^^^ exp { — dCQlIxjl^/fc}, 
thus for fc S {1 ,..., fco + f} we get 

p{x,k) < C||a:|r‘*exp{ - 2(70110:11^/(^0 + r)}|^ exp { - 2Go||a:||Vfc} 

< (7'||x||“‘^exp{ - (7o||a;||Vfco}- 
Hence, we can estimate 


/c=l 


Si{x,n) < C"\\x\\ ‘^expj - (7o||a;||^/fco} ^ 

For X £ and n G N, we define 
ko = 




4 log{nV’(||a:|| 2)}_ 

Then we have 

(5.4) 5i(a:,n) = o(n||a:||"V(lkir^)) • 

Claim 3. There are C > 0 and R > 0 such that for all n £ N and x £ Z'^, if 
||a:|| > R then 

(5.5) ko > C\\x\\'^f‘^, 

(5.6) nil}{kQ^) < (7(nV:(||x|r2))^^““^^^^^ 

Assuming for a moment the validity of Claim 3, we proceed with the proof of 
the theorem. We need to find the asymptotic of a sum 

r — l 

S2{x,n) = p(a:, fc + j)P(r„ = fc). 

k>ko j=0 

To do so we want to use the asymptotic from Corollary 2.1. Since (5.5) implies 
that for k > ko, 

Ikll < 

and 

k-^\\xf < Ck-'^l\ 

there is C > 0 such that 


y^p(x, fc + j) — r(27r) '^/^(det(5) ^^'^k exp{—||o:||^/(2fc)} 

3=0 


< (7fc-'^/2-2Aexp{-||a:||V(2fc)} 
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for all k > ko. Then, by setting 

I{x,n) = r(27r)“'^/^(det ^ exp{—||a;||^/(2fc)}P(T„ = k), 

fc>fco 

we obtain 

\S 2 {x,n) —I{x,n)\ < C ^ ||x||^/(2fc)}P(T„ = k) 

(5.7) < Cifep n). 

Next, for Gn{t) = 1 — Fn{t) = P(t„ > t) we can write 

(5.8) I(a;, n) = —r(27r)“‘^/^(det f exp{—||a:||^/(2t)} dG„(t). 

J [feo.oo) 

Let 

I{x,n) = r(27r)“‘^/^(det(3)"^/^ y ^ exp{-||a;||^/(2t)}^G„(t) dt. 

Since in (5.8) the integrand is continuous and G„ is monotone, the integration by 
parts yields 

|l(a:,n) - I{x,n)\ < exp{-||a;||^/(2fco)} 

(5.9) =o(||x|r‘^mA(||x|r2)), 

where in the last step we have used the definition of k^. To find the asymptotic 
behaviour of I(n,x) we compute 

(5.10) ±(t-^/^exp{-\\xf/i2t)}) 

= exp{-||x|| V(2t)} + exp{-||x|| V(2t)}. 

Next, let us define two integrals 

Ji = - / t"'^/^"^exp{-||xf/(2t)}G„(t) dt, 

J2 ^ 11^ f t-'i/2-2gxp|_||3^||2^^2i)}G„(t) dt. 

^ Jko 

Then, by (5.10), we have 

I{x,n) = r(27r)“‘^'^^(det(3)“^^^(T2 - Ji)- 
To show the asymptotic of Ji and J2 we need the following lemma. 

Lemma 5.4. For any /3 > 1 

^00 

/ t“^£(t)exp{-||a:||^/(2t)}G„(t) dt 

Jko 

^2/5W2-i£^+^,^^„||^||-2(/5-1)^(||^||-2^^ 
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as n'0(||x|| tends to zero. 

Proof of Lemma 5.4- First, we change variables by setting u = ||a;||“^t 


n||x||-2(^-i)^(||x||-2) Ao 


t ^exp{-||x||^/(2t)}G„(t) dt 

-0 GnjWxW^u) ^{\\X\\-'^U-^) 
l\\x\[-^ko nilj{\\x\\-^u-^) V-dkll"^) 

To estimate the integrand, by (5.6) and Theorem 2.2, we have 

Gn(lklpM) 


f 

J Ha; 


exp{—l/(2u)} du. 


n'ip[ ||a;||“2M“i) 


#(l |2 


<G, 


and, by (2.8) 


II /o r /3-1 Q-1 

' ’ < GM-“/ 2 j^ax{u-—,w—}, 


x||- 2 ) 


thus by applying the dominated convergence we obtain 

-B Gn{\\x\\^u) V'(lkll“^U“0 


x||-2fco 


nw\\\x\\~^u~ 


0 




exp{—1/(2'«)} dw 


pOO 

■ / M-'5-“/2exp{-l/(2M)}dM, 

Jo 


r(l - a/2) J, 

as n^(||x||“^) tends to zero, and the proof of the lemma is hnished. □ 

We now return to showing the asymptotic of I(x, n). For Ji, by applying Lemma 
5.4, we obtain 


(5.11) 


J, ^ ^ 2 ('^+“)/ 2 L[|^±^n||x||->(||x||- 2 ), 


as n^(||x|| tends to zero. In the case of J 2 , we obtain 


(5.12) 




r(l - a/2) 


as rM/;(||x|| ^) approaches zero. Putting (5.11) and (5.12) together we conclude 
(5.13) I{x,n) ~ rCd,an\\x\\~‘^'tp{\\x\\~'^). 

Now, by (5.7) and (5.9), 

S 2 {x,n) = {l{x,n) + o(n||x||“V(lkir^))) (l + o(l)), 
thus (5.13) yields 

S 2 {x,n) = n\\x\\~‘^'ip{\\x\\~'^) {rCd,a + o{l)), 

what together with (5.4) gives the desired asymptotic (5.3). 

To finish the proof of the theorem we need to justify Claim 3. First, by (2.8), 
there is i? > 0 such that for ||x|| > R and n G N 
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Hence, 

-log{nV'(||a:||"^)} < 3alog||a;||, 

and therefore 

fco > C'lkllVlogIkIh 

which easily implies (5.5). 

For the proof of (ii), again using (2.8), there are R> 0 and C > 1 such that for 
all n G N and ||a;|| > R 

K \ (l-a/2)/2 
-log{rM/;(||a;|r2)}j 

Therefore we obtain 

/„ / \ ( 1 + q / 2)/2 

^(^o) < CniAdla^lL )( - log {niAdlglL )}) 

< c (,.^dixr»))'‘-"™y 

and the proof of Claim 3 is finished. □ 

Remark 2. One possible application of Theorem 5.3 would be to study the Green 
function of the subordinate random walk 

G^{x) = ^p^{x,n). 

n>0 

We have 

G^{x) = ^p(x, k)G{k) 

k>l 

where {C{k) : fe > 0) satisfies 

k 

(7(0) = ! and G{k) =''^^c{tp,n)G{k — n), k>l. 

n—0 

Observe that the asymptotic formula of Theorem 5.3 does not provide enough 
information to obtain asymptotic formula for the Green function G^{x) at infinity. 
One way to study asymptotic decay of Gy, at infinity would be to study asymptotic 
properties of the sequence (G(fc) : fc > 0). Under certain special assumptions on 
^ this problem was solved in the recent paper of Bendikov and Cygan [3] . In this 
connection we also mention the closely related paper of Williamson [22]. 
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